A Back-Reaction Induced Lower Bound on the Tensor-to-Scalar Ratio 
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There are large classes of inflationary models, particularly popular in the context of string theory 
and brane world approaches to inflation, in which the ratio of linearized tensor to scalar metric 
fluctuations is very small. In such models, however, gravitational waves produced by scalar modes 
cannot be neglected. We derive the lower bound on the tensor-to-scalar ratio by considering the 
back-reaction of the scalar perturbations as a source of gravitational waves. These results show 
that no cosmological model that is compatible with a metric scalar amplitude of « 10~ 5 can have 
a ratio of the tensor to scalar power spectra less than « I0 -8 at recombination and that higher- 
order terms leads to logarithmic growth for r during radiation domination. Our lower bound also 
applies to non-inflationary models which produce an almost scale-invariant spectrum of coherent 
super-Hubble scale metric fluctuations. 

PACS numbers: 98.80.Cq 



I. INTRODUCTION 

Direct measurement of primordial gravitational waves 
or of the tensor to scalar ratio, r, would provide cosmol- 
ogists with invaluable information about the state of the 
inflationary universe - specifically, the value of the Hub- 
ble rate. The reason for this is that, unlike detecting the 
scalar amplitude, detecting the tensor amplitude unam- 
biguously determines the value of the inflationary scale 
since [l[ 
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where slow-roll conditions have been assumed. The ratio 
r is defined as the ratio of the power spectra of tensor to 
scalar modes. 

The current experimental upper bound is r < 0.5 and 
it is anticipated that, in the next few years, we will be 
able to probe for values of r in the range 1CP 1 — 10~ 2 Q. 
In the approximation in which both the scalar and tensor 
spectra are computed in linear cosmological perturbation 
theory, the value of r can be used to distinguish different 
classes of inflationary models. Whereas single scalar field 
models of large field inflation Q typically predict a value 
of r of around 10 -2 (plus or minus an order of magni- 
tude), small field models of inflation 0, [j| with poten- 
tials of Coleman- Weinberg type @ and hybrid inflation 
models Q typically predict a very small amplitude. For 
example, for a potential which in the slow-roll region can 
be approximated by 

V = V ~ A0 4 , (2) 
where A <C 1 is a self-coupling constant, the ratio is 

r~A 2 (^-) 2 . (3) 
mpi 

In the case of hybrid inflation, the small value of r can 
be seen either by direct computation, or else by making 



use of the "Lyth bound" 



r < 
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TV 2 mpi 
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where TV is the number of e-foldings of inflation between 
the time when cosmological scales exist the Hubble radius 
during inflation and the end of the inflationary phase, 
and A(jj is the change in the value of <p during the cor- 
responding period of inflation (which in the case of hy- 
brid inflation is much smaller than mpi. In particular, 
since the power spectrum of gravitational waves is given 
by the Hubble constant during inflation, low scale mod- 
els of inflation can have tiny values of r, values as low as 
10 -24 appearing in the literature [9(. Most of the recently 
popular string-motivated brane-antibrane inflation mod- 
els (see |l0|] for recent reviews) are of hybrid inflation 
type and hence lead to very small values of r, and in 
this context an upper bound on the primordial tensor to 
scalar ratio has been proposed [TTI | . 

The purpose of this short note is to point out that 
there exists an alternative mechanism for the production 
of tensors that is completely independent of the specific 
inflationary model. The back-reaction of scalar pertur- 
bations can act as a source of gravitational waves which, 
despite not being typically thought of as a main con- 
tributor to r in the case of large field inflationary mod- 
els, could dominate the value of r in other models. The 
main contribution to the induced gravitational waves is 
produced after reheating. Hence, our analysis and the 
resulting lower bound on r applies to all models which 
produce an almost scale-invariant spectrum of coherent 
super-Hubble scale fluctuations. 

Note that there has been a substantial body of previ- 
ous work on gravitational waves induced at higher order 
in perturbation theory in particular on gravi- 

tational waves produced at the end of inflation [3 [l5| . 
There has also been recent work on gravitational waves 
generated by the curvaton [16| . What is new in our work 
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is the focus on the modifications to the tensor to scalar 
ratio. Note that the conclusion that secondary gravita- 
tional waves may dominate over the primordial ones has 
also recently been reached in [l7l ]. 



II. GRAVITATIONAL WAVES FROM SCALARS 

Our starting point is the following perturbed line ele- 
ment written in terms of conformal time 77, and denoting 
the scale factor of the background by a{rj): 

ds 2 = a(r]) 2 l(l + 2<f>)drf-{l-2<S>)S l: jdx l dx : >+h lJ dx 1 dx j ]. 

The scalar metric fluctuations are described by a func- 
tion of space and time, and the gravitational waves are 
given by the transverse and traceless tensor hij. We are 
using longitudinal gauge to describe the linearized scalar 
metric perturbations [34[ . and are assuming the absence 
of anisotropic stress (see [19( for a comprehensive review 
of the theory of cosmological perturbations, and [201 f° r 
an introductory overview). For each wavenumber k, there 
are two polarization states for gravitational waves, each 
described by an amplitude h and a polarization tensor 
e ij ■ 

In other words, each second rank tensor hij can be 
decomposed as 

hij (r),x) = h (1) (r), x)e^ + h {2) (77, x)ef j } . (6) 

In what follows, we make use of this decomposition to 
express the dynamics of the gravitational waves, hij, in 
terms of the scalar coefficient functions h^' 2 '. 

At linear order in the fluctuations, the scalar and ten- 
sor modes are independent. However, if the magnitude of 
the linearized gravitational waves is very small compared 
to the magnitude of the scalar modes, then the secondary 
tensor fluctuations generated by non-linear interactions 
from the scalar modes cannot be neglected. We can de- 
termine the coupling between $ and hij by expanding 
the gravitational action in powers of the perturbation 
amplitude. At leading order, we find: 

Sgrav = ^ [ d 4 X^R (7) 
1D7T J 

= ^Jd 4 x[^f((hr-(dm 

-a 2 (r/)/i(8$V 2 $ + 6 (<%$) 2 ) + ...] . 

In the above (and below), we have supressed the indices 
denoting the graviton polarization state as both states 
couple identically to matter. 

Additional sources for gravitons arise from the matter 
sector coupled through the kinetic term. In the spirit 
of simple inflationary universe models, we will take the 
matter sector to be described by a single canonically nor- 



malized scalar field <j). In this case we write 

Smatter = J d*X^ y/^gW d^tf - V (if)] (8) 

= J d 4 x... -ia 2 (r7)M(W 2 
+2<$>d i (pd l 6ip + 2<f> 2 (d l tp) 2 ) ... , 

where we have extracted only the terms which are linear 
in h. If the background is homogeneous, then the terms 
proportional to diip are absent. 

Combining ([7]) and ([5]) leads to the following, sourced 
equation of motion for the tensors: 

m" + (-V 2 - (9) 
a 

= -a(r/)[16$V 2 $ + 12(di$) 2 

O 

+ — (diSif 2 , 
m Pl 

where we have rescaled the tensor mode, h — ^i/a(vj). 
This transformation has the advantage of eliminating the 
Hubble damping term and making the particle produc- 
tion process more transparent. 

In the standard case, squeezing of the quantum vacuum 
state leads to particle production [2l| with the squeezing 
being determined by the dynamics of the background. 
The squeezing plays the role of a tachyonic instability, 
which can be seen by the presence of a negative mass 
term in ([§]). Exactly the same thing happens in the scalar 
sector except that, instead of the squeezing term being 
proportional to — — , it is proportional to — — with z = 

cijt- It is this different dependence on the behaviour of 
the inflaton that makes it possible to have a scenario in 
which the tensor amplitude is far below that of the scalar. 

We now consider a scenario in which the squeezing 
produces a negligible amount of gravitons, but a large 
amplitude of scalar metric fluctuations. In this case, we 
can determine the amplitude of the (secondary) tensors 
being produced via interactions from the scalar modes (in 
the language of particle physics, one could say "by scalar 
decay"). To this end, we neglect the squeezing term in 
© entirely. Since we focus on super-Hubble scale fluc- 
tuations, we can also neglect the spatial gradient terms 
[35| . With these approximations, the partial differential 
equation ([9]) reduces to an ordinary one which takes the 
form 

fi" = ia(ry)[16$V 2 $ + 12(5i$) 2 

+^(^M 2 ]- (10) 

The above equation ([1^) makes it clear that scalar met- 
ric fluctuations source gravitational waves at all times. 
We can simplify the equation by noting that the scalar 
sector has only one independent degree of freedom . 
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During slow-roll inflation we can make use of the con- 
straint equation [22j 



to recast (llOl) as 



5<p * -2^*, 



-a(r])[4$V 2 $ 



(11) 



(12) 
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After reheating, a formula analogous to (fTT))) but with 
matter taken to be a perfect fluid must be used. 

We first focus on the gravitational waves generated 
during inflation. Making use of the definition of the slow- 
roll parameter 



16tt 1 V ' 



(13) 



in (|12p . we see something perhaps a little surprising: the 
second term in the source (the term originating from the 
matter sector) scales as 1/e. As mentioned above, in 
the standard case (ignoring back-reaction) , a small value 
of e suppresses the amplitude of the tensors. However, 
including the effects of back-reaction, we see that a small 
e leads to a higher amplitude of induced gravitational 
waves. 

We now wish to use the above formalism to compute 
the amplitude of the secondary gravitational waves at 
the time when the mode of interest is crossing the Hub- 
ble radius during the initial inflationary phase, assuming 
that the initial state of all fluctuation modes on ultravi- 
olet scales is the quantum vacuum state. Because of our 
assumption on the initial state, it is only scalar modes 
which have already exited the Hubble radius and thus 
undergone some squeezing which can produce such sec- 
ondary gravitational waves. 

One way to compute the source term in (| 1 2[) is to ex- 
pand $ in Fourier modes, insert into the expression on 
the right hand side of the equation, use the fact that 
the scalar spectrum is scale-invariant on super-Hubble 
scales, and to extract the fc'th Fourier mode of /j, by in- 
verse Fourier transform of the result and by contracting 
with the polarization tensor of the gravitational wave. 
The integrals over momenta run from an infrared cut- 
off (which can be taken to be the scale corresponding 
to the Hubble radius at the beginning of inflation, but 
which does not play any role in this calculation) to mo- 
menta corresponding to the Hubble radius. Modes on 
sub-Hubble scales are in their vacuum state and hence 
cannot generate any gravitational waves. Also, as shown 
in [l5| . the effect of short wavelength scalar fluctuations 
would average to zero in our coupling term. 

A short cut to the calculation is to simply to promote 
$ to the status of quantum operator [lj| and to take the 
vacuum expectation value of the right hand side of (fT2l) , 



making use of the two-point function (in the limit that 
r -> 0) 



\®(x, r))<b(x + r, 7?)|0) 



dk sin(fcr) , .. l9 
T - fcr ' *' ' 



(14) 



and restricting the integral to the IR phase space, k = 
[0,7i], where 7i is the Hubble rate in conformal time. 
In the above, |<5$| 2 is the power spectrum of the scalar 
modes (including the phase space factor k 3 ). Multiplying 
with the polarization tensor to extract the amplitude /i 
of the gravitational wave, we thus get 
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through which the evolution of /i is now directly related 
to the scalar power spectrum |i5$| 2 . Note that in the 
above formula, the scalar metric power spectrum is the 
one during inflation. 

Several remarks need to be made about this calcula- 
tion: to obtain this result, we assumed scale- invariance 
of the scalars. As well, f| 15[) can only be used to extract 
information about /i on large scales, i.e. H 3> k p h ys - 
The physical reason for this condition is obvious - mo- 
mentum conservation forces the produced tensors to have 
smaller momenta than the scalars from which they orig- 
inate. Considering contributions arising solely from the 
scalar IR sector is not as limiting as it may seem due to 
the fact, through the course of inflation, the IR phase 
space grows exponentially while the UV phase space re- 
mains constant (see the arguments presented in |22j | and 
[23l]). However, in our calculation there is no IR diver- 
gence and the non-gradient contributions of the IR modes 
(those which in the case of adiabatic matter fluctuations 
are pure gauge 0, [H, [26[) do not contribute. Note 
that the - dependence on the size of the effects of back- 
reaction has also been noticed in [l2j , [13] , and in a very 
different context .[28|. 

Integrating (|15| immediately yields 



a(rj) 



1, 



(16) 



(17) 



where we have used a(rj) = — and 7i = a(rj)H. The 
scale invariance of the above is manifest. 

The amplitude of the power spectrum of scalar metric 
perturbations during inflation is suppressed compared to 
the spectrum after inflation by the factor e 2 , the inverse of 
the squeezing factor of $ 2 during the transition between 
the inflationary phase and the radiation phase of Stan- 
dard Cosmology. Taking the current power spectrum to 
be w 10~ 10 , we estimate the contribution to the tensor- 
scalar ratio induced by back-reaction before reheating to 
be 
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(18) 
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This is negligible compared to the contribution to r of 
the linear gravitational waves produced from quantum 
vacuum fluctuations. 

Next, we consider the gravitational waves induced by 
back-reaction after reheating. To obtain a lower bound 
on the effects of back-reaction, we neglect the matter 
contribution on the right hand side of (fT0|) . Hence, there 
is no e _1 enhancement factor. However, the amplitude 
of <!> to be used is not suppressed by e. Hence, repeating 
the steps which led to (fTT|) leads to an analogous formula 
without the e _1 term in the square parentheses on the 
right-hand side of the equation. This give a contribution 

r w 1CT 10 . (19) 

Hence, the back-reaction of scalar metric fluctuations has 
a more important effect after reheating than before. 

Our analysis is not restricted to the inflationary era. 
Interestingly, (|15[) . predicts the creation of gravitons on 
super-Hubble scales even after reheating. Direct substi- 
tution of the radiation-domination scale factor (a(rf) = 
77/770) leads to the result 

h( v ) ~ ln(a(77)). (20) 

At the onset of matter domination (a(r>) ~ f] 2 ), the 
amplitude of the super-Hubble tensor modes again be- 
comes constant. 

Taking the above into consideration, we take as our 
final estimate for the back-reaction determined tensor- 
to-scalar ratio at the time of recombination to be 

r a 1(T 8 . (21) 

At the order in perturbation theory we are now work- 
ing in we must also consider the decay of the gravitational 
waves into scalars through interaction terms of the form 

S h 2^ oc J d A xa 2 (^h] V 2 hj . (22) 

However, for the values of r obtained here, this decay is 
negligible. 

III. CONCLUSION 

We have estimated the amplitude of gravitational 
waves produced via the back-reaction of scalar metric 



fluctuations, focusing on their effect on the tensor-to- 
scalar ratio r. In conclusion, our result (fT9f shows that 
back-reaction of scalar modes will produce a sizeable 
gravitational wave background in a wide variety of in- 
flationary models in which the value of r computed from 
purely linear theory is very small. 

The contribution to r from back-reaction of scalar 
modes picks up its major contribution after reheating. 
Hence, our lower bound on r will be valid in any theory 
producing a roughly scale-invariant spectrum of fluctua- 
tions at late times, provided these fluctuations are coher- 
ent outside the Hubble radius (which appears to be de- 
manded by the data based on the acoustic peak structure 
of the observed angular power spectrum of cosmic mi- 
crowave anisotropies) . Examples of such theories are the 
Ekpyrotic universe scenario [29| , the Pre-Big-Bang model 
|3fJ | and the recently proposed string gas structure forma- 
tion paradigm [3lj . Our analysis is particularly relevant 
for the Ekpyrotic scenario in which a negligible amount 
of gravitational waves are generated in the linear theory. 
Back-reaction, however, will induce a contribution which 
is similar in magnitude to what is induced in many small 
field inflationary models. 

An interesting point that was noted was the ability of 
backreaction to increase the amplitude of IR modes after 
inflation. This is in stark contrast with the linear result 
in which the amplitudes of both scalar and tensor modes 
are frozen until the time of second Hubble crossing. 

This work is complementary to an earlier study of 
the back reaction effects of scalar modes onto the scalar 
modes themselves [32j which were shown to be negligi- 
ble. As this note was being finalized, a paper appeared 
[33l ] which studies vector modes produced by primordial 
density fluctuations. 
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